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Abstract 

We reconsider the Volume Source Technique (VST) for the determination 
of flavor singlet quantities on the lattice. We point out a difficulty arising in 
the case of fermions in real representations of the gauge group and propose 
an improved version of the method (IVST) based on random gauge transfor- 
mations of the background configuration. We compare the performance of 
IVST with the method based on stochastic estimators (SET). We consider the 
case of the N=l Supersymmetric Yang-Mills Theory (SYM), where just one 
fermionic flavor is present, the gluino in the adjoint representation, and only 
flavor singlet states are possible. The work is part of an inclusive analysis of 
the spectrum of the lightest particles of the theory, based on the simulation of 
the model on a 16 3 • 32 lattice with dynamical gluinos in the Wilson scheme. 
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1 Introduction 



Supersymmetry (SUSY) is broken on the lattice owing to the finite lattice spacing a. We 
consider the N=l Supersymmetric Yang-Mills Theory (SYM) with gauge group SU(2) 
and Wilson discretization in the fermion sector. Here SUSY is also explicitly broken by 
the Wilson term. However, by properly tuning the (renormalized) gluino mass to zero, 
SUSY is expected to be recovered in the continuum limit pQ with exponentially small 
O(a) deviations. 

Manifestation of SUSY occurs at the non-perturbative level, the most interesting phe- 
nomenological implication being the expected ordering of the bound-states of the theory 
in supermultiplets. In the low energy sector in particular, effective Lagrangians for SYM 
predict |2J E] two Wess-Zumino supermultiplets. The spin-0 particles are represented 
by meson-like bound states of the gluino and by glueballs, respectively, of opposite par- 
ity (this classification is of course only valid in absence of mixings, which are however 
expected). The spin- t> particle of the multiplet is in both gluino-glue bound-state. 

We focus here on the problem of determining the masses of meson-like gluino bound 
states. Borrowing the terminology of QCD, these represent "flavor singlet" states. Indeed, 
SYM resembles Nf=l QCD, with the quark in the fundamental representation replaced 
by the gluino in the adjoint representation. The lattice computation of flavor singlet cor- 
relators is difficult because of the presence of disconnected diagrams (see [3] for a recent 
review on the topic). The exact evaluation of the correlator for these diagrams is not 
feasible since it requires the trace over color and space-time indices of the fermion prop- 
agator in the background of the gauge configuration, which in turn involves the solution 
of an "all-points to all-points" inversion problem for any given gauge configuration. The 
first approach to the subject was based on a volume source jHj, the so-called "Volume 
Source Technique" (VST). For a given background configuration the method delivers an 
estimate of the correlator which, however, contains spurious terms represented by non- 
closed loops. In j^j, where QCD was considered, it was argued that these terms disappear 
in the ensemble-average on the basis of gauge invariance. In this paper we reconsider this 
argument more generally, showing that it is not applicable to models where the fermions 
are in real representations of the gauge group, as is the case for any representation of 
SU(2) and for the adjoint representation of SU(iV c ). We propose a new formulation of the 
method, based on random gauge transformations of the background gauge configuration, 
which solves the problem. We shall refer to this in the following as to IVST (Improved Vol- 
ume Source Technique). Due to the randomness introduced by the gauge transformation, 
IVST is analogous to the well known Stochastic Estimator Technique SET [6 . In both 
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cases the systematic error introduced by the computational procedure is converted into 
a statistical one and can be controlled by increasing the number of stochastic estimates. 
As a consequence IVST and SET can be directly compared. 

This work represents the sequel of a long-standing project having the goal of a lattice 
verification of the non-perturbative low-energy properties of SYM. We refer to [Zj and the 
references therein for the scope and goals of past studies. The model is simulated by means 
of the dynamical-gluino two-step multi-bosonic algorithm. Details on the algorithm can 
be found in [H]. The present analysis is based on a sample of configurations of SU(2) SYM 
on a 16 3 • 32 lattice. Partial results have been reported in [Oj. 

In the next section we shall reconsider the theory of VST and propose the improved 
version of it, IVST. In section 3 the numerical results will be presented, comparing IVST 
and SET; finally section 4 contains our conclusions. 



2 The Volume Source Technique revisited 

In this section we consider lattice gauge theory with gauge group SU(iV c ). The results of 
primary interest are for the gauge group SU(2) or for models with fermions in the adjoint 
representation of the gauge group. This includes SYM in particular. In the following, 
Greek letters denote Dirac indices, Latin letters color, Tr^ and Tr c are the respective traces. 
With the usual bilinears il){x)Til)(x) as insertion operators for the singlet mesonic states, 
where T = 1 or 75, the disconnected part of the mesonic correlator in the background of 
a gauge configuration {U} can be written as 

C r , d isc[U}(x - y ) = ^-Tr d [TS(x )] Tr d [TS(y )} , (1) 

where the time-slice sum S(xq) represents the trace over color and space indices of the 
inverse fermion-matrix, i.e. the propagator in the background of the gauge configuration 

{uy. 

S aP (x ) = J2^c[Q^J. (2) 

X 

VST delivers an estimate of S a p(x Q ) at the price of a single inversion for each value of the 
color and Dirac index. The inversion problem with the volume-source uy reads 

qz = cV a,al , (4 a,al w = s ab <W (3) 

with solution 

^xb/3 = [Q U v' \xbj3 — Q x b/3,xaa + ^xb/3,yaa • (4) 

y^x 
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When Z^- a '°^ in the above equation is used to estimate the time-slice sum (J2J), 



(5) 



the last term in Q yields contributions to the disconnected part of the correlator (0) 
which represent non-closed loops. Such elements of the inverse fermion-matrix with x ^ y 
are non-gauge-invariant and are canceled in the average over the gauge-ensemble (which 
is gauge- invariant). However, there are also contact terms in the correlator, which are 
potential sources of systematic errors. 

In the original work [S], which introduced VST in the context of QCD, these unwanted 
terms were avoided by considering the correlator 



with one of the time slices conjugated. Owing to the fact that the product 3 <8> 3 of 
fundamental representations of SU(3) does not contain the trivial representation, a gauge 
invariant contact term does not appear. The argument holds more generally for the 
fundamental representation of SU(iV c ) for N c > 2. 

In the case of gauge group SU(2), which has real representations only, or in the case of 
the adjoint representation of SU(iV c ), this prescription, however, does not help. For SU(2) 
the product of two fundamental representations contains the trivial one, which leads to 
non-vanishing contact terms again. The same is true for the adjoint representations of 



We now want to consider the gauge invariance of the contact terms in detail. We focus 
on the correlator for © the discussion is analogous. 

Consider the following average over gauge transformations g(x) (gauge-average): 



CrMUKxo - yo) = ^Tr d [rS(x )} Tr d [r^(y )] 



(6) 



SU(AQ. 





The gauge- average induces an average over the gauge-orbit {U 9 }. Using 



QxlP 9 ] = g\x)Q- x l[U}g{y) 




and the general formula 




fundamental 



adjoint 
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(in the adjoint representation g are real orthogonal matrices of dimension N% — 1), the 
gauge-average of Q reads for xo ^ yo 

(S a p(x )S jS (yo)) 9 = Tr c [Qi J ^ Tr c [Q~^ yS ] +AJ2 ^c[Q^, y pQ^J ■ (10) 

x y x,y 

The above expression represents the gauge-invariant part of S a p(xo)S^s(yo)- 

Let us now consider the ensemble-average of S a p(xo)S^s(yo)- In the limit of infinite 
statistics any given gauge-orbit is completely covered, implying that the ensemble-average 
delivers in particular a gauge-average. Using the result in (JTUJ) this implies 

(s a(3 (x )S^(yo)) u = (S a p{x )S 7 s{yo))u + A \J2Tr c [Q-^ yl3 Q~^ xS } J . (11) 

\ x,y I v 

We thus obtain that replacement (0) in (JTJ) produces an error-term for the full disconnected 
correlator 

Cr,diBc(% - yo) = C r ,di SC (x - y ) + ACr idi8C (a; Q - Vo) , (12) 
AC r , disc (a;o - yo) = Al/^Tr c [Tr d [g-Jr]Tr d [g^r]]\ . (13) 

S \ x,y I u 

The conclusion is that the error-term in (j3J) produces a systematic error in the correlator, 
which does not vanish in the ensemble- average even in the limit of infinite statistics. This 
error is due to gauge invariant contact terms in the correlator, as shown above. The 
spurious term resembles the connected contribution 



C v ^ nn [U}{x -y,) = ~fyY^Ti cd [Q-lTQ y lT}, f=\ 

S x.y I ' 



fundamental , , 
(14) 

adjoint 

;v v 
the only difference being in the Dirac structure and the numerical factor. This outcome 
is not surprising considering that gauge invariance strongly constrains the space-time and 
color structure. We have checked the presence of the error-term numerically for both 
types of correlators ((TJ) and (jBJ) for gauge group SU(2), see sectional 

At this point we make the simple observation that the error is removed by using the 
gauge-average of S a p(x ) to determine the time-slice sums, since 

(^S a /3(xo)j = S a p(x ) . (15) 

In practice this is obtained by averaging S a p(xo) over a sufficiently large number N g of 

gauge configurations obtained from the original one by random gauge transformations [H] 1 

1 After the completion of this study we noticed that the use of random gauge transformations in VST 
was recently pointed out in 0]. 



g(x), namely with a flat probability distribution 

4-1. (16) 
dg 

where dg denotes the Haar measure on the gauge group. Besides solving the problem of 
the error (JT3j) in the correlator, the method brings the additional benefit of disentangling 
the systematic error inherent in VST from the statistical one: in the limit of an infinite 
number of random gauge transformations N g — > oo the former goes to zero, only the 
second one surviving. In this view the improved version of VST (IVST in the following) 
is analogous to the techniques based on stochastic estimators (SET), the randomness of 
the source being replaced by that of the gauge transformation. 2 This allows for a direct 
comparison of the two methods, which is carried out in the next section. 



3 Numerical analysis 

The simulation parameters of the gauge sample are (3 = 2.3 and n = 0.194. The estimated 
value of the lattice spacing is, in QCD units, a « 0.06 fm (a -1 « 3.3 GeV); there are 
indications JU] that the gluino is still relatively heavy (m§ > 200 MeV on the basis of 
QCD-inspired arguments). The set-up of the two-step multi-bosonic algorithm algorithm 
is the same as in [TJ], and ~4000 thermalized configurations were stored every 5 or 10 
cycles. In order to obtain an estimate of the autocorrelation time of the disconnected 
part of the mesonic correlator, an analysis of the autocorrelation time of the smallest 
eigenvalue of the hermitian fermion-matrix was performed. The procedure is based on 
the expectation that the disconnected part of the mesonic correlator is strongly related to 
the infrared behavior of the fermion-matrix. After that, a subsample of 218 supposedly 
uncorrected configurations was selected. This constitutes the sample for the numerical 
analysis. 



3.1 Time-slice sums 

For each configuration, 50 estimates of the time-slice sums (j2J) were performed, each ob- 
tained by applying a random gauge transformation on the original gauge configuration as 
explained in the previous section. The computations were performed in 64 bit arithmetic. 
Improved summation techniques were employed to ensure accuracy. 

2 Actually on the basis of © IVST could be seen as a stochastic estimator method with a particular 
stochastic volume source. 
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In the case of SYM the Majorana nature of the gluino field (invariance under charge 
conjugation) allows to compute the inverse of the fermion- matrix for only half of the 
matrix-elements in Dirac space. This implies that, in the case of SU(2) SYM, only 6 
fermion-matrix inversions must be performed for each configuration, compared to 12 inver- 
sions needed for QCD. So the total number of inversions N{ nv required for a determination 
of the time-slice sum with N est estimates is N{ nv = 6N est . 3 

As IVST is based on stochastic estimations, a comparison with stochastic-source meth- 
ods SET suggests itself. We consider the SET variant with complex Z 2 noise in the spin 
explicit variant SEM |12j . In this case each estimate of the time-slice sum is obtained 
by inverting the fermion-matrix with source (wg )xbp — dapVxb where ry^ are indepen- 
dent stochastic variables chosen at random among 775 (±1 ± i)- For SET one has then 
N inv = 2N est . (Again a factor of two less comes from the symmetry of SYM.) We computed 
165 estimates of the time-slice sums, in this case using 32-bit arithmetic. 

In Fig. [T]the evolution of the estimated value of Tr[Q _1 r] = J2 Xo r ^ r d[S(xo)T] for a 
chosen configuration is displayed as a function of the number of needed inversions N{ nv . 
The error bounds represent the statistical uncertainty on the stochastic estimation. For 
both IVST and SET the value stabilizes after 150-200 inversions, with compatible results. 
This test on a single configuration only serves as a cross-check of the two methods, the 
physical information being contained in the ensemble-averages, Fig. 121 In the scalar case 
the two methods give compatible results after only 50 inversions. In the pseudoscalar case, 
fluctuations much larger than the error-bounds indicate additional effects. The fluctua- 
tions appear to be more relevant for SET, where 32-bit arithmetic was used. Moreover, 
in the latter case the estimate has an offset, while in the case of IVST the expected value 
(zero) is approached after ~100 inversions. 

The evolution of the statistical error of the estimation for one configuration is displayed 
in Fig. El showing the a priori non-obvious result that the two methods introduce the same 
amount of stochastic uncertainty. The error in the estimation of the ensemble-average is 
shown in Fig. HJ We see that in both cases the error stabilizes after 100 inversions. In 
the pseudoscalar case, IVST seems to out-perform SET, although the large instabilities 
prevent us from drawing firm conclusions. 

3 N es t coincides with N g of previous section. The change of notation is for the sake of the homogeneity 
when comparing with SET. 
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3.2 Correlators and masses 

In order to show the effect of the error-term (|13|). we computed the disconnected correlator 
in two ways: i) following the correct procedure according to Eq. (|T3|) (IVST); ii) perform- 
ing the gauge-average as in (JTUJ) . As one can see in Fig. El for the pseudoscalar meson, the 
error-term produces a sizeable effect on the disconnected correlator. IVST and SET are 
in good agreement. The effective mass is shown in Fig. The impact of the error on the 
effective mass is suppressed in the first time-slices where the connected contribution (JTlj) 
plays a larger role. However in the last time-slices, where the disconnected contribution 
dominates, the effect of the error-term shows-up in the form of a pronounced instability 
of the effective mass as a function of the time-separation (for At=13 an estimate is not 
even possible). In the last few time-separations At = 14, 15, IVST delivers a better result 
compared to SET (no estimate is possible with SET for At = 15). Since the disconnected 
contribution to the mesonic correlator is essentially of infrared nature, the region of large 
time-separations is important for the determination of masses. 

4 Conclusions 

We propose an improved version of the Volume Source Technique (IVST) which elimi- 
nates erroneous contact terms in the case of fermions in real representations of the gauge 
group. The improved version is based on random gauge transformations and is analo- 
gous to stochastic estimator methods (SET). Comparison between IVST and SET shows 
agreement and substantial equivalence. In few cases, e.g. for the determination of effec- 
tive masses, IVST seems to give slightly better results. A study with higher statistical 
precision should put these observations on firmer ground. 
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Figure 1: Evolution of the estimated value of Tr[Q _1 ] and Tr[Q _1 75] for a chosen config- 
uration as a function of the number of the needed inversions (with error-bounds). Full 
lines: IVST, dashed lines: SET. 
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Figure 2: Evolution of the average value of TrfQ" 1 ] and Tr[<5 _1 75] over the complete 
sample as a function of the number of the needed inversions (with error-bounds). Full 
lines: IVST, dashed lines: SET. 
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Figure 3: Evolution of the statistical error of the estimated value of TrfQ -1 ] and Tr[Q _1 7 5 ] 
for the same configuration as in Fig. 1, as a function of the number of the needed inversions. 
Full lines: IVST, dashed lines: SET. 
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Figure 4: Evolution of the statistical error on the average value of Tr[Q _1 ] and Tr[Q _1 75] 
over the complete sample as a function of the number of the needed inversions. Full lines: 
IVST, dashed lines: SET. 
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Figure 5: The disconnected pseudoscalar correlator Cr,disc(At). 




Figure 6: The effective mass of the pseudoscalar meson. 
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